Band structure of ZnO from resonant x-ray emission spectroscopy 
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Soft x-ray emission and absorption spectroscopy of the O K-edge are employed to investigate 
the electronic structure of wurtzite ZnO (0001). A quasiparticle band structure calculated within 
the GW approximation agrees well with the data, most notably with the energetic location of the 
Zn 3d - O 2p hybridized state and the anisotropy of the absorption spectra. Dispersion in the band 
structure is mapped using the coherent k-selective part of the resonant x-ray emission spectra. We 
show that a more extensive mapping of the bands is possible in the case of crystalline anisotropy 
such as that found in ZnO. 

PACS numbers: 71.55.Gs, 78. 70. En, 78. 70. Dm, 71. 15. Mb 



I. INTRODUCTION 

The wide band gap semiconductor zinc oxide (ZnO) 
is a material with much potential in future electronic 
devices. 1 Despite significant interest in its fundamental 
electronic properties^ few spectroscopic studies of the 
full band structure exist. A small number of angle- 
resolved photoemission spectroscopy (ARPES) studies 
have been reported ^2*2*4 yielding some agreement with 
density functional theory-based (DFT) calculations and 
providing information on surface and defect states. How- 
ever, ARPES is limited to probing occupied states, and 
charging effects, surface preparation and sample qual- 
ity are serious issues that can limit the accuracy of re- 
sults. For example, photoemission measurements have 
variously located the fully occupied Zn 3d semi-core lev- 
els at energies ranging from 7.5 to 8.8 eV below the Fermi 
level jSAMaSaSi The close proximity of the 3d level to the 
O 2p-derived valence band has an appreciable impact 
on band structure calculations f£ & 1Q i n i 12 i 13 As a result, 
there is a need for bulk sensitive measurements of the 
ZnO shallow core-level and valence band dispersion. 

Here, we present angular dependent resonant x-ray 
emission spectroscopy (RXES; also known as resonant 
inelastic x-ray scattering - RIXS) of the O K-edge of 
wurtzite ZnO crystals, for comparison with hybrid DFT 
calculations and quasiparticle energies obtained within 
the GW framework. Earlier studies used x-ray ab- 
sorption and emission spectroscopy (XAS and XES) 
to directly measure the O 2p partial density of states 
(PDOS) of nano-structured ZnQ 14 i 15 and its anisotropy 
in the conduction band. 16 We go further by exploiting 
the dipole- and k-selection rules of RXES to map the 
anisotropic valence band of single-crystal ZnO. 

RXES band mapping has been reported in simple ma- 
terials like graphite, BN, SiC, and more j 17 i 18 i 19 and re- 
ports on more complicated compounds related to ZnO 



are appearing j 20 i 21 i 22 i 23 However, to our knowledge there 
is no literature examining the K-edge electronic struc- 
ture of a post transition metal oxide in the detail we 
present, or on any material using our image-based tech- 
nique. Further, we show crystal anisotropy to be a useful 
tool, allowing a larger part of the Brillouin zone to be 
uniquely accessed via RXES. Finally, compelling exper- 
imental evidence for use of the final-state rule in RXES 
is presented. 



II. EXPERIMENT AND THEORY 

The sample, a 500 nm ZnO epilayer, was grown on 
epi-ready (0001) sapphire substrates by plasma- assisted 
MBE, as described previously. 24 As the sample was part 
of a separate study on doping, approximately 0.01% Ag 
was incorporated during growth, resulting in a layer av- 
eraged electron concentration of 7.3 x 10 16 cm -3 . The 
high crystalline quality of the film was confirmed by RBS 
under the channeling condition^ a streaky RHEED pat- 
tern, and by low temperature photoluminescence, which 
shows donor bound exciton peaks are dominant with up 
to three LO phonon replicas. 

The x-ray spectroscopy was performed on the undu- 
lator beamline X1B at the National Synchrotron Light 
Source at Brookhaven National Laboratory, which is 
equipped with a spherical grating monochromator and 
a Nor dgren- type emission spectrometer. The energy res- 
olution over the O K-edge was approximately 0.20 eV 
and 0.37 eV for the XAS and RXES respectively. XAS 
were recorded in total electron yield mode. The photon 
energy was calibrated using the O K-edge and Ti L-edges 
of rutile Ti02 measured during the experiment. The XES 
was calibrated with the 2nd order L-edge of Zn from both 
a calibration metal sample and the ZnO sample. 

Due to its wurtzite structure, the ZnO p-projected 
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FIG. 1: (color online) p xy (black) and p z (red) orbital con- 
tributions to XES (closed symbols) and XAS (open symbols) 
compared to PDOS from the GW calculation (lines). Bars 
indicate RXES excitation energies. 



PDOS is strongly anisotropic^ To measure the 
anisotropy XAS and RXES were recorded for light in- 
cident at 20° and 70° relative to the sample normal, re- 
ferred to as normal and grazing geometry, respectively. 
Dipole selection rules mean that XAS recorded in the 
normal (grazing) geometry is dominated by the contri- 
butions from O 2p xy (2p z ) orbitals (see inset, Fig. [I]). 
The opposite is true for RXES as the emission spectrom- 
eter is oriented perpendicular to the incident light: x-rays 
incident normal (grazing) to the sample result in grazing 
(normal) RXES. Thusp xy RXES is analysed withp^ XAS 
and vice versa. This correspondence has been observed 
for wurtzite GaN by Strocov et al.. 20 

Electronic structure calculations were performed 
within the hybrid DFT (HDFT) framework using the 
HSE03 functional for exchange and correlation. 26 Quasi- 
particle effects were taken into account by a subsequent 
GW correction of the HSE03 eigenvalues using many 
body perturbation theory. In the GW calculations, the 
Coulomb potential was fully screened using the random 
phase approximation dielectric function based on the 
HSE03 eigenvalues and functions^ It has been shown 
that the ab initio combination of HDFT and GW calcula- 
tions gives excellent values for the fundamental gaps and 
(i-band positions of many semiconductors. 28 . At present 
calculating a fully k-resolved GW band structure is too 
computationally intensive, so we compare XAS and XES 
to PDOS obtained from the GW calculation, and RXES 
to a modified HDFT band structure. 



III. RESULTS AND DISCUSSION 

Figure [I] shows the O K-edge XAS (open symbols) 
and XES, taken well above threshold (closed symbols). 



FIG. 2: (color online) Comparison between GW (red line) 
and HDFT (solid black line) calculated O p-projected PDOS 
of the occupied states. The GW calculation more accurately 
positions the Zn 3d-0 2p hybrid when compared with exper- 
iment (closed symbols). Multiplying the HDFT energy scale 
by a factor of 1.1 results in improved agreement (dashed line). 



Strong linear dichroism associated with the anisotropy 
in the unoccupied states is seen in the XAS, especially 
in the first 10 eV above the conduction band minimum 
(CBM) where the p xy (black) contribution is consider- 
ably stronger than the p z (red) contribution. The XES 
is comprised of O 2p states near the valence band max- 
imum (VBM) and a clear peak corresponding to O 2p 
states hybridized with Zn 3d orbitals at 519.9 eV. 

These spectra are compared with PDOS (solid lines) 
obtained from the GW calculation. The PDOS have 
been convolved using gaussians with widths of 0.37 eV 
and 0.20 eV for the occupied and unoccupied states re- 
spectively. This is an estimate of the intrinsic instrumen- 
tal resolution of the experiment; no lifetime effects have 
been considered. The PDOS have been aligned so that 
the main valence band peak matches with experiment. 
This yields excellent agreement with the entire XES, ac- 
curately locating the 3d hybrid. The weak anisotropy 
in the PDOS, combined with the lower resolution of the 
XES, explains the lack of contrast between p xy and p z 
XES. The final state in XAS includes a core hole that is 
not included in the calculation, and which prevents ab- 
solute comparison between the energy scales of the XAS 
and XES*22 Therefore, to compare to XAS, the theoreti- 
cal unoccupied PDOS has been rigidly shifted by an ad- 
ditional — 1.0 eV, providing an estimate of the core hole 
binding energy. With this adjustment the experimental 
agreement is excellent in terms of both the main peak 
locations and intensities, and the orbital anisotropy. 

The success of the GW calculation is highlighted by 
a comparison with the HDFT result. We calculate a 
band gap of 3.2 eV, compared to 2.1 eV obtained from 
the HDFT calculation (and 3.4 eV experimentally^) The 
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FIG. 3: (color online) HDFT band structure. Calculation is 
orbitally resolved into p xy (black) and p z (red) contributions. 
Error bars represent strength of contribution. The energy 
axis has been scaled for better agreement with GW theory 
(see Fig. [2}. Labels indicate RXES excitation energies. 
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FIG. 4: (color online) Normalized coherent RXES showing 
p xy (black) and p z (red) data. Labels indicate excitation en- 
ergy (hi/). 



valence band PDOS is also greatly improved: Figure [2] 
shows the total O 2p-projected occupied PDOS calcu- 
lated using HDFT (black line), and with the GW correc- 
tion applied (red line). With GW the Zn 3d - O 2p hy- 
brid peak is shifted to lower energy, peaking at — 7.1 eV, 
in much closer agreement with the experimental peak at 
—7.4 eV (closed symbols), correcting a key flaw of stan- 
dard DFT calculations. 28 With the Zn 3d states located 
lower in energy the rest of the valence band is widened 
towards lower energies due to reduced p-d repulsion. We 
note that scaling the HDFT energy axis by a factor of 1.1 
is enough to obtain relatively close agreement with the 
GW calculation (dashed line) for the occupied states. An 
additional upward shift of 0.5 eV to the conduction band 
is sufficient to achieve similar agreement for the unoc- 
cupied states. We have confirmed the accuracy of this 
approach approach by calculating the GW band struc- 
ture at selected k points. 

Figure [3] shows the band structure obtained from the 
HDFT calculation, with the above corrections applied. 
The p xy (black) and p z (red) orbital contributions to the 
eigenvalues are plotted independently; error bars repre- 
sent the relative weight of the contribution. To simplify 
RXES analysis, the p xy (p z ) unoccupied states are plot- 
ted with p z (p X y) occupied states (see above). RXES can 
be viewed as a scattering process where the final state 
includes an electron in the conduction band and a hole 
in the valence band, but no core-hole. So, unlike Fig. [IJ 
no core hole correction is applied to the conduction band. 

As expected, the band structure is that of a direct-gap 
semiconductor. The strongly dispersing band that forms 
the CBM at T is also anisotropic: only p xy states are 
found along T-M and T-K; p z along T-A-H. This combi- 



nation of strong dispersion and orbital selectivity make 
ZnO especially suitable for band mapping with RXES. 
For weakly correlated systems, and excitation energies 
close to the absorption edge, resonant effects contribute 
to XES, adding a k-conserving coherent contribution to 
the scattering. The crystal anisotropy, combined with 
orbital selection rules, adds further k-selectivity to the 
RXES. For ZnO, the use of these properties in concert 
allows valence band dispersion to be measured across L- 
M-r-A-H-K-r of the Brillouin zone. 

RXES were taken over a range of excitation energies, 
indicated by vertical bars on Fig. [TJ and horizontal bars 
on Fig. 03 The coherent fraction, extracted using stan- 
dard techniques, 19 ranged from 0.45 to 0.6. 

Using the above we construct Figure O which shows 
normalized intensity maps of the coherent fraction of each 
spectrum placed in k-space according to the position of 
the bars of Fig. [3j Figures [SJa) and (b) show the p xy and 
p z data, respectively. The imaging approach highlights 
the information contained in the raw spectra. Dispersion 
toward higher energy at T, as expected for a direct-gap 
material^ is immediately obvious for both states. For 
the p xy states dispersion along T-A-H-K is clearly re- 
solved, with the upper band dispersing to lower energy. 
The bands around —4 eV decrease in energy and increase 
in intensity along A-H-K, in line with theory. The posi- 
tion of the Zn 3d states is apparent, just below the theo- 
retical bands. The spectrum obtained with an excitation 
energy of 537.6 eV is located near both M and K, out 
of order with the other spectra, because of the strongly 
localised states in the conduction band at 10.5 eV. The 
p z states add M-T and K-T to the parts of the brillioun 
zone we can resolve. The lower energy band near —4 eV 
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FIG. 5: (color online) The coherent RXES overlaid with or- 
bitally resolved band structure: (a) p xy and (b) p z orbitals. 
Intensity scales from zero (white) to one (black). 

disperses to —5.5 eV approaching T, in direct contrast 
to the observed behaviour of the p xy band in this energy 
range, and in agreement with theory. 

There has been much discussion about the correct 
use of intermediate and final states in RXES calcula- 
tions.— i2L2£i We note that if the core- hole correction were 



applied to the unoccupied states of Fig. [3] (a downward 
shift of the theory by 1.0 eV) the p xy spectrum taken with 
537.6 eV would be relocated to near T and the strong in- 
tensity around —4.5 eV would be anomalous. This is 
clear evidence for the use of the final-state rule in RXES. 



IV. CONCLUSION 

In conclusion we have used the strong anisotropy and 
dispersion of wurtzite ZnO, the k-selectivity of the RXES 
technique, and an HDFT+GW calculation to construct a 
band mapping of ZnO across a wide range of high symme- 
try points. The energy of the Zn 3d core level is located 
and evidence for the use of the final state rule in RXES 
is presented. It should be possible to use further high 
resolution RXES experiments to tune theoretical calcula- 
tions across the technologically important post transition 
metal oxide and nitride series. 
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